Collective dynamics of complex plasma bilayers 
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A classical dusty plasma experiment was performed using two different dust grain sizes to form 
a strongly coupled asymmetric bilayer (two closely spaced interacting monolayers) of two species of 
charged dust particles. The observation and analysis of the thermally excited particle oscillations 
revealed the collective mode structure and wave dispersion in this system; in particular the existence 
of the theoretically predicted k — energy (frequency) gap was verified. Equilibrium molecular 
dynamics simulations were performed to emulate the experiment, assuming Yukawa type inter- 
particle interaction. The simulations and analytic calculations based both on lattice summation and 
on the QLCA approach are in good agreement with the experimental findings and help identifying 
and characterizing the observed phenomena. 
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Particle bilayers (parallel planes occupied by interact- 
ing particles and separated by a distance comparable to 
the intcrparticle distance within the layers) can be viewed 
as an intermediate stage between two-dimensional (2D) 
and three-dimensional (3D) systems. It is the interplay of 
the 3D interaction and the 2D dynamics that creates the 
rich new physics predicted and observed in interacting bi- 
layers that makes these systems interesting in their own 
right. At the same time, bilayer configurations are also 
ubiquitous in widely different physical systems. Of spe- 
cial importance are those of charged particles (with like 
charges: unipolar bilayer, or with opposite charges: bipo- 
lar bilayer). Examples are semiconductor heterostruc- 
tures [if, cryogenic traps 0], overdamped system of lipid 
membranes [J, interfacial superconductors [if, etc. 

From the theoretical point of view, during the past 
three decades unipolar layered systems were studied in 
the weak coupling limit by means of analytic calculations 
[B|, [f| and in the strongly coupled regime by semi-analytic 
lattice calculations |7|, approximate liquid state calcula- 
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tions 8| and computer simulations [9j, [lfj, lll|, 
More recently, bipolar (electron-hole) bilayers have been 
the focus of intense computer simulation efforts, both 
in the classical 14 , lB|, EH] and in the quantum regimes 
171 1 X 81 ] , where features such as structural phase transi- 
tions, bound dipole formation [l||, etc. were detected. 

Very recently a seminal observation by Hyde et. al. [2(| 
has led to the realization that strongly coupled unipolar 
bilayers can be created in laboratory complex (dusty) 
plasma environments. This can be accomplished by us- 
ing a mixture of two differently sized grains: in view 
of their necessarily different charge to mass (Z/m) ra- 
tios, the two species would settle at different equilibrium 
heights in the plasma sheath, as governed by the local 
balance of gravitational and electric forces. Thus this 



novel type of bilayer would be, in contrast to most of the 
previously observed ones, a binary bilayer with hitherto 
unexplored features. Its structural properties were al- 
ready reported in [2(| ■ Subsequent computer simulations 
predicting collective excitations (wave propagation) and 
their dispersion were carried out in recent years both for 
Coulomb and Yukawa type isotropic interactions [2l], tHJ . 

In this Letter we report on the experimental inves- 
tigation of the collective dynamical properties of a bi- 
nary dusty plasma system. Our results constitute the 
first observation of the mode spectrum of a strongly cou- 
pled (liquid or solid) bilayer. Earlier results |23j were 
restricted to the weakly coupled state, where all collec- 
tive modes exhibit an acoustic behavior. We confirm the 
predicted [Tj, |8[ benchmark of the strongly coupled mode 
structure, the development of optic modes with a wave- 
number k — "energy (frequency) gap" . 

Our dusty plasma experiments have been carried out in 
a custom designed vacuum chamber with an inner diam- 
eter of 25 cm and height of 18 cm. The lower, powered, 
18 cm diameter, flat, horizontal, stainless steel electrode 
faces the upper, ring shaped, grounded aluminum elec- 
trode with an inner diameter of 15 cm at a height of 
13 cm. Experiments have been performed in (4.6 purity) 
argon gas discharge at a pressure p — 0.8 ± 0.05 Pa, in 
a steady gas flow of ~ 0.01 seem, with 13.56 MHz radio 
frequency excitation of ~ 5 W power. In the experi- 
ment melamine-formaldehyde micro-spheres with diame- 
ters d\ = 3.63±0.06/jm and d,2 = 4.38±0.06/jm are used. 
For illumination we apply a 200 mW, 532 nm dpss-laser. 
Our CCD camera has a resolution of 1.4 Megapixels and 
runs at 29.54 frames per second acquisition rate. Particle 
masses are: m\ = 3.8 10 -14 kg and 777.2 = 6.6 10 -14 kg. 

During the evaluation of the raw images (typically over 
60000 per experiment) identification and position mea- 
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surement of t he p articles is performed using the method 
described in [24|. The identification of the light and 
heavy (smaller and larger) particles is based on their 
scattered intensities, which is in this size domain pro- 
portional to the square of the diameter. After tracing 
the particles motion from frame to frame we obtain the 
positions and velocities of each particle as a function of 
time. The layer separation was measured by simply ro- 
tating the camera and illuminating setup and taking im- 
ages through the side window. After calibrating with 
a size standard we find the following average structural 
parameters: number of particles: N\ — 682, N2 = 636; 
the Wigner-Seitz radius calculated from the total num- 
ber of the observed particles (N) and the area (A) of the 
field of view: a = y/A/wN = 0.243 mm (in the follow- 
ing distances appear normalized to a); layer separation 
d = d/a = 0.43. 

Based on the particle positions the g(r) pair distri- 
bution function is obtained and compared to molecular- 
dynamics (MD) simulation results, (which was carried 
out by using the same input parameters as in the exper- 
iment), as shown in Fig. [TJ Comparing experiment and 
simulation, one can conclude, that both capture qualita- 
tively identical features: order, shape and position of the 
peaks agree satisfactorily. Differences in the amplitude 
and decay rate are due to friction and density gradients 
due to the finite and confined nature of the experimen- 
tal system. In the simulation infinite system size (peri- 
odic boundaries) and no friction are assumed. Also, the 
layer assignment procedure during the data processing 
has some degree of uncertainty. Peak positions are con- 
sistent with the underlying hexagonal lattice structure, 
in agreement with theoretical results in 0, [H| • 

Visual observation of the recorded images confirms 
that most of the system is in a hexagonal configuration, 
sometimes small domains with rhombic and square unit 
cells were also found. Although the ground state config- 
uration for a bilayer is a regular staggered rectangular 
lattice, a large degree of substitutional disorder [25j , due 
to the small inter-layer separation and the finite temper- 
ature, was observed, as shown in the particle snapshots 
in Fig. ED 

In the MD simulation Yukawa type pair-interaction 
§ij = (e 2 Z i Z J /4:'Ke )eyi-p[-ri j /X D ]/r i: j, where r i3 is the 
three-dimensional distance between particles i and j, is 
assumed using ~ 4000 particles with periodic boundary 
conditions. Since the interlayer separation is of the or- 
der of 1/4 Ac, the attractive force between grains in the 
two layers due to ion focusing or wakefield effects are as- 
sumed to be small, since these effects occur generally at 
larger separations (see e.g. [lij]). The neglect of such 
attractive forces is further bolstered by the experimen- 
tal observation that the grains in the two layers do not 
form strings, but rather assume a staggered structure. 
Particle densities, masses and layer separation are taken 
from the experiment, k = a/Xjj = 0.5 and Z\ = 2550, 
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FIG. 1: Pair-distribution functions and particle snapshots 
obtained from the experiment (a,c) and the MD simulation 
(b,d). Indexes and symbol types label the layers [open: 
layer 1, filled: layer 2]. 



Z-2 = [dajd\)Z\ » 3080 arc further assumed; this param- 
eter set was found in lattice calculations to best repro- 
duce experimental dispersions. 

In the experiment we observe a pressure jump of about 
38% when switching on the discharge. Taking into ac- 
count the active and the total volume of the vacuum 
system, this results in an « 85% average temperature 
increase in the discharge region, resulting in T p « 550A'. 
Measuring the dust particle velocities gives the direct ki- 
netic energy, resulting in T v i « 400-ftT and T V 2 ~ 5QQK 
for the two layers, and the Coulomb coupling parameters 
[defined as T m = (e 2 /^e ){Z 2 Jk B T m a m )]: T 1 » 800 
and ri ~ 900. To obtain dynamical information on 
the system's collective excitations we use the method 
based on the Fourier transform of the microscopic density 
and current fluctuations, already successfully applied in 
MD simulations, 



From the Fourier trans- 



see e.g. [11|, l27|. Knowing the par- 
ticle positions vs. time, first we calculate the micro- 
scopic densities, as well as longitudinal and transverse 
currents: p (m) (k,t) = J2j exp(-ik • rj), A< m >(k,i) = 
Ej v j,\\ exp(-ik • Tj), T (m )(k,t) = J2j V j-X exp( 
for the two layers (to = 1,2 
forms of p(k, t) — * p(k,o>), A — > A(k, w), and r — > r(k, u>) 
one can calculate the power spectra, e.g. S^nfk, ui) oc 
(p^i— k, — u>)p( n >(k, uj)), averaging is over the time- 
slices available. In the present case of an asymmetric 
bilayer the labeling of the modes is not possible in a sim- 
ple way, as it is for symmetric bilayers llj, where the 
two "+" and "-" polarizations clearly separate and cor- 
respond to in-, and out-of-phase oscillations. Therefore, 
in the following we do not label the modes, only indicate 
in which spectra they appear as peaks. For an analysis of 
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the mode structure one has to examine the spectra indi- 
vidually and identify the peak positions. Sample spectra 
are shown in Fig.[5]to illustrate the fundamental features 
in more detail. Plotted are results of the bilayer experi- 
ment and of the MD simulation; also to serve as reference 
standard, results of a separately performed single layer 
experiment carried out in the same experimental setup 
and under the same conditions except for using only one 
particle species (with diameter 4.38 ± 0.06/im, the to- 
tal particle number in the field-of-view was N — 1945, 
resulting in higher density, thus higher nominal plasma 
frequency as in the bilayer case). 
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FIG. 2: Sample spectra illustrating the principal features of 
(a) Ln longitudinal, (b) Tu transverse current fluctuations 
and (c) Si 2 inter-layer density fluctuations. The one-layer 
spectra from the experiment and MD simulations are shown 
and compared to experimental single-layer spectra obtained 
in a separate measurement. The inter-layer spectra (exper- 
iment and MD simulation) show a positive and a negative 
peak representing the in-phase and out-of phase oscillations, 
respectively. The vertical scale is shifted for clarity. 

In Fig. HJa) and (b) one can observe both in the ex- 
periment and in the simulation the presence of a strong 
primary peak at lower frequencies and a shoulder-like, 
weak and wide feature at higher (lu « 50s" 1 ) frequen- 
cies. The strength of this high-frequency peak becomes 
more obvious in comparison with the single layer situa- 
tion, where at the same frequency the spectral power has 
already dropped two orders of magnitude, while it has 
the same slope as seen in the bilayer spectrum at higher 
frequencies (u > 70s -1 ). Figure [2jc) shows the inter- 
layer density fluctuation spectrum S12 at ka — 0.4 from 
the experiment and the MD simulation. The novel fea- 
ture here is the appearance of a negative peak indicating 
an out-of-phase oscillation at the higher frequency, coin- 
ciding with the "shoulder" in the Lu and Tu spectra. 



As it is expected, all spectral peaks are more pronounced 
and sharper in the simulation, where frictional damping 
and disturbing effect of the finite confinement are absent. 

We compare our observational results with the mode 
dispersion calculated through two theoretical models. 
The first model is a perfect staggered rectangular lattice, 
with phonon propagation along the two principal axes. 
The lattice calculation is based on the formalism used 
in [3] for electronic bilayers, adopted to the asymmetric 
Yukawa system. The second model is a completely disor- 
dered solid bilayer, with the mode dispersion calculated 
in the Quasi-Localized Charge Approximation (QLCA) 
formalism. The QLCA formalism was adapted to the 
binary Yukawa bilayer system (based on [23, H^l) and 
the dispersions for the 4 modes were calculated with the 
input of the experimental pair-correlation functions. It 
should be noted that while for symmetric configurations 
the (in general) 4x4 dynamical matrix can be decomposed 
into two 2x2 matrices, resulting in a clear separation of 
the collective modes into longitudinal vs. transverse and 
in-phase vs. out-of-phase polarizations, in the present 
case, this kind of mode separation, in general, is not pos- 
sible. It can be done only along the directions of the prin- 
cipal axes of the crystal, or in an approximation where 
the system is assumed to be isotropic as in the QLCA 
description. 




FIG. 3: (color online) Lu one-layer longitudinal (a,c) and Tu 
one-layer transverse (b,d) current fluctuation spectra. Peak 
positions (hot colors) mark the dispersion co(k). Black sym- 
bols in (a,b) are frequencies form lattice summation calcula- 
tions including the two principal lattice directions. Lines in 
(c,d) are the corresponding QLCA dispersions. 

Figure [3] shows an example of the longitudinal and 
transverse current fluctuation spectra in the upper layer, 
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color coded in the wave-number (ka) j frequency (to) 
plane. Overlayed are a selected set of lattice dispersions 
(panels [a] and [b]), and the QLCA mode dispersions 
(panels [c] and [d]). An additional, very low frequency, 
linear dispersion can be seen, which is most likely the fin- 
gerprint of an unavoidable net motion of the ensemble. 



30 - + 



20 ,J / 



la) 



2 3 
ka 



2 3 
ka 



.-;%H, 



FIG. 4: (color online) Symbols: peaks identified in the lon- 
gitudinal (a) and transverse (b) experimental spectra. Lines: 
selected lattice dispersions with x and y polarizations. 

Figure [4] shows the more detailed dispersion properties 
of the dusty plasma bilayer together with lattice calcu- 
lation results. The modes are not labeled, for reasons 
discussed above, only their sources (the specrta in which 
the peak was found) are indicated. 

In Figs. [3] and Hone can observe the expected acoustic 
modes similar to those in 2D Yukawa layers, and an opti- 
cal excitation with finite and nearly constant frequency, 
even at very low wave-numbers. This latter mode is iden- 
tified as the so called energy (frequency) gap, studied in 
great detail in earlier simulations and theoretically pre- 
dicted through lattice calculations and for strongly 
coupled liquids through the QLCA approach [§, 11|. The 
QLCA model provides the value of the gap frequency as 



gap 



Z2112 Z 2 mi 
2 V Z\ n\ Z\m,i 



T(r)g 12 (r)Ar, (1) 



where uif — (e 1 Z\n\) j {2e^m\^ a\a 2 \ the kernel T(r) — 
z- 3 re- w [(l + w)(l - 3d 2 /z 2 ) + w 2 ], z = yjr 2 + d 2 , and 
w = z/Xd. 

In view of the prevailing local lattice structure in co- 
existence with a high degree of disorder it is not a pri- 
ori clear, which of the two models (lattice vs. QLCA) 
should provide a better description of the mode struc- 
ture. Inspection of Fig. [3]shows that for low k values the 
acoustic portions of the low frequency modes are equally 
well described by either model. For higher k values the 
repeated Brillouin zone structure is clearly visible, indi- 
cating the superiority of the lattice model. For the high 
frequency optic mode the QLCA predicts a single "fre- 
quency gap" at k = with Eq.©: = 41.4s" 1 and 
w gap = 43.5s" 1 calculated with the input of experimen- 
tal and simulation gi2(r) data. This value and the high-fc 



tapering off of the optic modes seem to be more along the 
line of the QLCA description. 

We can conclude, that our dusty plasma experiment 
using two different dust sizes created a strongly coupled 
bilayer system that can be well approximated by a unipo- 
lar binary Yukawa bilayer model. This model has served 
as the basis for MD simulations, lattice calculations, and 
a QLCA calculation. All these approaches are in good 
agreement with the experiment, and verify the presence 
of an optical collective mode characterized by a finite en- 
ergy (frequency) gap at k = wave-number, distinguish- 
ing the strongly coupled bilayer from a weakly coupled 
one, where all the modes have an acoustic character. 
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